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ABSTRACT 

Transits on single stars are rare. The probability rarely exceeds a few per cent. Further¬ 
more, this probability rapidly approaches zero at increasing orbital period. Therefore transit 
surveys have been predominantly limited to the inner parts of exoplanetary systems. Here we 
demonstrate how circumbinary planets allow us to beat these unfavourable odds. By incorpo¬ 
rating the geometry and the three-body dynamics of circumbinary systems, we analytically 
derive the probability of transitability, a configuration where the binary and planet orbits 
overlap on the sky. We later show that this is equivalent to the transit probability, but at an 
unspecified point in time. This probability, at its minimum, is always higher than for single 
star cases. In addition, it is an increasing function with mutual inclination. By applying our 
analytical development to eclipsing binaries, we deduce that transits are highly probable, and 
in some case guaranteed. For example, a circumbinary planet revolving at 1 AU around a 
0.3 AU eclipsing binary is certain to eventually transit - a 100% probability - if its mutual 
inclination is greater than 0.6°. We show that the transit probability is generally only a weak 
function of the planet’s orbital period; circumbinary planets may be used as practical tools 
for probing the outer regions of exoplanetary systems to search for and detect warm to cold 
transiting planets. 

Key words: binaries: close, eclipsing, spectroscopic - astrometry and celestial mechanics: 
celestial mechanics, eclipses - planets and satellites: detection, dynamical evolution and sta¬ 
bility, fundamental parameters - methods: analytical, numerical, statistical 


1 INTRODUCTION 

In the burgeoning search for extra-solar planets, circumbinary plan¬ 
ets represent some of the most exotic systems found to date. They 
pose astronomers with interesting questions regarding their de¬ 
tectability (Schneider 1994), abundance (Armstrong et al. 2014; 
Martin & Triaud 2014), formation (Pierens & Nelson 2013; Kley 
& Haghighipour 2014), habitability (Haghighipour & Kaltenegger 
2013; Mason et al. 2014), orbital dynamics (Leung & Hoi Lee 
2013) and stability (Dvorak 1986; Dvorak et al. 1989; Holman & 
Wiegert 1999). 

Answers to these questions are reliant on planet detections. So 
far there have been reported discoveries from several techniques, 
including PSR B1620-26 (Thorsett et al. 1999 with pulsar timing), 
HD 202206 (Correia et al. 2005 with radial velocimetry), DP Leo- 
nis (Qian et al. 2010 with eclipse timing variations), Ross 458 (Bur- 
gasser et al. 2010 with direct imaging), and Kepler-16 (Doyle et 
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al. 2011 with transit photometry). There are presently ten transit¬ 
ing circumbinary planets known, all found by the Kepler telescope 
(Welsh etal. 2014). 

The advantage of finding circumbinary planets in transit is that 
they can yield an unambiguous detection, thanks to a unique sig¬ 
nature that is hard to mimic with false positives. The photometric 
measurement of the radius can be complemented with transit tim¬ 
ing variations (TTVs), eclipse timing variations (ETVs) or spec¬ 
troscopy to obtain the mass and bulk density, which are important 
from a formation perspective. Transits also open the door to atmo¬ 
spheric characterisation (Seager & Deming 2010), the measure of 
the Rossiter-McLaughlin effect (Queloz et al. 2000; Fabrycky & 
Winn 2014), and the detection of exomoons (Kipping et al. 2012). 

It will be shown in this paper that circumbinary planets, be¬ 
yond their exoticity, are useful astronomical tools. Their particu¬ 
lar geometry and orbital dynamics lead to potentially much higher 
transit probabilities in comparison with single stars. There is also a 
weaker dependence on orbital period, allowing us to extend transit 
studies to the outer regions of stellar systems. 

The paper is structured as follows. In Sect. 2 we introduce the 
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geometry of circumbinary planets. Next in Sect. 3 we analyse the 
orbital dynamics of circumbinary systems and the effects on their 
observability. We then define the concept of transitability and ana¬ 
lytically derive a criterion for its occurrence in Sect. 4. Following 
this, we convert this criterion into the probability of a circumbinary 
system exhibiting transitability in Sect. 5, similar to the work done 
for single stars (Borucki & Summers 1984; Barnes 2007). In Sect. 6 
we analyse the special case of eclipsing binaries. 

As an observer, the observable quantity is a transit, not tran¬ 
sitability. This is why in Sect. 7 we connect the two concepts, ver¬ 
ifying that a system exhibiting transitability is effectively guaran¬ 
teed to transit, albeit at an unspecified point in time. Some illus¬ 
trative transit wait times are calculated, revealing that they may be 
within a few years for many systems. In Sect. 8 we discuss some ap¬ 
plications and limitations of our work, before concluding in Sect. 9. 


2 GEOMETRY 

We will treat a circumbinary system as a pair of Keplerian orbits 
in Jacobi coordinates, with the addition of first-order dynamical ef¬ 
fects (Sect. 3). The inner orbit is the stellar binary (subscript “bin”). 
The outer orbit is the planet around the binary centre of mass (sub¬ 
script “p”). Each Keplerian orbit is an ellipse characterised by four 
orbital elements: the semi-major axis a, eccentricity e, argument 
of periapsis tu and true anomaly /. These quantities are defined in 
Fig. 1. This set of four is not unique, and often we will use the pe¬ 
riod T instead of the semi-major axis. The orientation of each orbit 
in three dimensions is defined using two extra angles: the inclina¬ 
tion 1 and longitude of the ascending node fl. In Fig. 2 we depict the 
3D orientation of a binary and planet orbit. We take the observer to 
be looking down the z-axis. An eclipsing binary, for example, cor¬ 
responds to /bin ~ 7r/2. Throughout this paper we use radians unless 
otherwise specified with a ° symbol. 

The orientation of the planetary orbit with respect to the binary 
is characterised by two quantities: the mutual inclination 

cos AI = sin /bin sin /p cos Afl + cos /bin cos /p, (1) 

and the mutual longitude of the ascending node 

An = nbin - np, (2) 

which are also shown in Fig. 2. When using transit photometry or 
radial velocimetry, the observer is sensitive to AQ but not to the 
individual quantities nbin and np. Throughout this paper can we 
therefore take nbin = 0 and allow np to vary. 


3 DYNAMIC ORBITS 
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Figure 2. A circumbinary planet in a misaligned orbit (blue, outer) around 
a binary star' system (pink, inner). The misalignment is characterised by the 
mutual inclination. A/, and the mutual longitude of the ascending node, AH. 
The observer is looking down the z-axis from above, and hence the grey x-y 
plane denotes the plane of the sky. 


the binary plane (AD = 0 to 2;r), whilst maintaining A/ = const 
(Schneider 1994; Farago & Laskar 2010; Doolin & Blundell 201 1). 
The precession period Tprec according to Schneider (1994) is 


A static Keplerian orbit is insufficient for accurately describing a 
circumbinary planet. Owing to perturbations from the binary, the 
orbital elements defined in Sect. 2 vary on observationally rele¬ 
vant timescales. We include in our derivation the most prominent 
of these effects: a precession in the planet’s orbital plane. This be¬ 
haviour is described by a time-variation in A/, and Af2. We restrict 
ourselves to circular binaries and planets*. In this case, the orbital 
plane of the planet rotates at a constant rate around the normal to 




1 

cos AI ’ 


(3) 


where the stars are assumed to be of equal mass. An alternative, 
more complex derivation can be found in Farago & Laskar (2010). 
The planet orbit is stable as long as it is not too close to the binary. 
An approximate criterion from the work of Dvorak (1986); Dvorak 
et al. (1989); Holman & Wiegert (1999) is 


' It is technically a misnomer to speak of circular circumbinary orbits, since 
perturbations from the binary cause gp to vary even if initially zero. How¬ 


ever, this only has a very small effect on the transit geometry. See Sect. 8.5.2 
for further detail. 
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Figure 3. Surfaces of section of the mutual inclination. A/, and mutual lon¬ 
gitude of the ascending node, Af2, between the binary and planet orbital 
planes. 


flp ^ Sflbm- (4) 

In Fig. 3 we demonstrate how precession effects A/ and AO 
using numerical N-body integrations". We ran a set of simulations 
with Af2 starting at 90° and A/ varied between 0° and 180° in steps 
of 10° and, each corresponding to a different curve in Fig. 3. The 
stars are of mass \Mq and Q.SMq with abin = 0.07 AU. The planet 
is a massless test particle with = 0.3 AU. The green, inner curves 
are for prograde orbits with clockwise precession. The blue, outer 
curves are for retrograde orbits with anti-clockwise precession. The 
gap between the green and blue curves corresponds to A/ = 90°, 
i.e. for a polar orbit where the precession period becomes infinitely 
long (Eq. 3). 

These orbital dynamics have observational consequences. The 
inclination planet on the sky, /p varies with time according to 

/p = A/cosj;^rj-r/bin, (5) 

\ prec / 

where t is time. An example is shown in Fig. 4 for two circumbi¬ 
nary systems. The binary in both systems has equal mass stars with 
Ma = Mb = IMo, flbin = 0.1 AU, h\n = 110°, and Dbm = 0°. 
The planet is a massless body with starting values /p = 130° and 
flp = 0°. The mutual inclination is 20°. The two planets shown in 
the figure have different values for Op) 0.3 AU for the black, dashed 
sinusoid and 0.6 AU for the blue, solid sinusoid. 

The maximum and minimum values of /p are independent of 
flp. The planet semi-major axis does, however, strongly influence 
the precession period: Tp^c = 5.6 yr for = 0.3 AU and 65.9 yr 
for flp = 0.6 AU, according to the N-body simulation. The analytic 

^ All N-body simulations in this paper are done using a fourth-order 
Runge-Kutta algorithm, where energy loss due to its non-symplectic nature 
was kept to negligible levels. 



Figure 4. Variation of /p over time for two circumbinary systems (Op = 0.3 
AU in black and white dashes, Op = 0.6 AU in blue). The horizontal red 
line denotes the binary’s orbital plane inclination on the sky, /bin- The grey 
region corresponds to when the planet is in transitability. 



expressions in Eq. 3 produces precession periods of 6.0 yr and 67.1 
yr, showing it to be reasonably accurate. 


4 CRITERION FOR TRANSITABILITY 

Transitability is an orbital configuration where the planet and bi¬ 
nary orbits intersect on the sky, like the example shown in Fig. 5. 
In this scenario transits are possible but not guaranteed on every 
passage of the planet past the binary, because of the relative motion 
of the three bodies. This terminology was first introduced in Martin 
& Triaud (2014), we where formally defined and elaborated upon 
a concept that had been already used in several studies (Schnei¬ 
der 1994; Welsh et al. 2012; Kratter & Shannon 2013). In Fig. 4 
the grey region denotes the time spent in transitability. During each 
precession period, there will be zero, one or two intervals of tran¬ 
sitability, or permanent transitability in only two scenarios: 1) /bi„ 
and /p are both very close to n/2 and 2) polar orbits where /yn = 0 
and /p = n/2. 

We work to derive a criterion that predicts whether or not a 
circumbinary planet will enter transitability at any point during the 
precession period. As a first approximation, we know that the planet 
is in transitability when the planet orbit is perpendicular to the plane 
of the sky: 
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l,."y (6) 

This is the most conservative case possible, since it ignores the fi¬ 
nite extent of the binary. According to Eq. 5, /p is guaranteed to 
reach njl if 


A/ > 


n 

2 



(7) 


This is the first-order criterion for transitability. Three things are 
apparent: 1) mutual inclinations aid transitability, which is contrary 
to conventional views on transit geometries, 2) this criterion is in¬ 
dependent of the planet period and 3) this criterion is easiest to ful¬ 
fil at /bin ~ tt/2, i.e. for eclipsing binaries. This criterion was also 
derived by Schneider (1994), who was the first author to analyse 
circumbinary transit probabilities in the presence of precession. 

The second level of complexity is to include the full extent of 
the stellar orbits, meaning that a value of /p offset from njl may 
still exhibit transitability. Consider the limiting case of transitabil¬ 
ity. This is when the planet and binary orbits barely overlap when 
|/p -;r/2| is at a minimum {dl^ldt = 0). We calculate the orientation 
of the binary and planet orbits in this configuration. Take Eq. 1 and 
rearrange it to isolate the term containing Afl: 


cosAf2 : 


cos A/ - cos /bin cos /j 
sin/bin sin/p 


( 8 ) 


In Eq. 8 only /p and AH are time-dependent quantities. Differenti¬ 
ating both sides of Eq. 8 with respect to time leads to 


- sin AD- 


.dAQ. 

dt 


dip 

cos /bin sm /p — sm l^n sm /p 
dL 


- (cos A/ - cos /bin cos L ) sin /^n cos /p — 

dt 

/(sin^ /bin sin^ 7p)- 


(9) 


By substituting dL/dt = 0 into Eq. 9 we get 


sinAD—;—= 0. (10) 

dt 

From Sect. 3 it is known that dASlIdt = const t 0, and hence 
sin AD = 0, implying that AD = 0. This means that the limiting 
case of transitability occurs when the ascending nodes of the binary 
and planet orbits are aligned. This simplifies the geometry and cal¬ 
culations. According to Sect. 3, AD is guaranteed to equal zero at 
some point during the precession period. When AD = 0 the mutual 
inclination calculation in Eq. 1 is simplified to 


A/ = |/M„-/p|. (11) 

In Fig. 6 we show a circumbinary system in the limiting case 
of transitability. This is a “side-on” view of the orbit, with the ob¬ 
server located to the right of the page. Using this diagram, we define 
half the projected heights on the sky of the two stellar orbits to be 


T/a.b — flA.B sin 



+ Ra.b , 


( 12 ) 


where 


(13) 



Figure 6. A side-on view of a circumbinary system in the limiting case of 
transitability, where AD = 0. The two extreme vertical positions of each 
star on the sky are drawn in different colours. In this example, the planet is 
barely in transitability on the secondary star, but not on the primary. 


are the individual semi-major axes for the two stars and 
_ A1b,a 

/^b.a - -r. - — 

Ma -b Mg 

are the reduced masses. Similarly for the planet. 


Xp = flp sin 



There is transitability on stars A and/or B when 


(14) 


(15) 


2fp<XA,B, (16) 

According to the orbital dynamics (Eq. 5), the planet will enter tran¬ 
sitability at some point (fulfilling Eq. 16) if the following criterion 
is met: 


AI > 


n 

2 


i^A-B (7bin)i 


(17) 


where the angle y0A,B is a function of the binary extent on the sky, as 
seen by the planet. This is the second-order criterion for transitabil¬ 
ity. Combine Eqs. 11 and 17 to get 


l^bin lp\ 



n 

1 > 

2 ~ 


/iA.B(7bin)- 


(18) 


By inserting Eqs. 12 and 15 into Eq. 16 and rearranging to match 
the form of Eq. 18, we obtain 


y8A,B(7bi.i) = sin 


<3a,B 

71 

^bin 

Ra.B 

_l. _ 

dp 

2 

dp 


(19) 


where the quantities inside the square brackets are sufficiently small 
that we can use the small angle approximation. 

The size of jSa.b determines how easy it is for a system to ful¬ 
fil the transitability criterion in Eq. 17. Depending on the separate 
values of /3a and /3b, it is possible to fulfil the criterion for just one 
of the stars. Generally /3b > /3a except for eclipsing binaries. 

To test the validity of the second-order transitability criterion 
we ran numerical N-body simulations on thousands of hypothetical 
circumbinary systems. The details are shown in Appendix A. The 
analytic criterion is shown to be very accurate, with an error less 
than 0.1%. All error cases were near the limit of the inequality in 
Eq. 17. Errors arise due to small variations in the semi-major axis 
and eccentricity, which in the limiting case of transitability may 


Oa.B - <tbiii/^B,A 
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lead to a contrary result to the prediction of Eq. 17. This is elabo¬ 
rated upon in Sect. 8.5.2. 


5 PROBABILITY OF TRANSITABILITY 

For a given set of orbital parameters, we can calculate the proba¬ 
bility that a given observer will observe transitability, at some time 
during the precession period. As an initial approximation, we use 
the first-order transitability criterion (Eq. 7). The orientation of a 
circumbinary system on the sky is uniformly random. It therefore 
follows that cos /yn has a uniform distribution, and hence the prob¬ 
ability density function is 


P(fbin) “ Stn/bin. (20) 

By integrating this between the bounds specified by Eq. 7, we ob¬ 
tain an approximate probability of transitability: 

f A.B “ I sin /bin^fbin 

= sin(A/). (21) 

The probability is period-independent and non-zero, except for 
strictly coplanar systems. 

The next step is to include the finite extent of the binary or¬ 
bit. The probability of transitability is a function of the size of the 
solid angle subtended on the celestial sphere such that the planet 
and stellar orbits are seen overlapping, including the full orbital 
evolution. We demonstrate this in Fig. 7a. To simplify the diagram, 
we only draw the angles for transitability on the primary star, but 
the calculation proceeds identically for the secondary. The planet 
orbit is shown in two different positions, corresponding to the two 
extrema of /p separated by 2A/ (Eq. 5). The angles ^a.b and 0 a.b 
are the angles subtended by the planet in transitability on each star. 
They are functions of how big the stellar orbit is, as seen by the 
planet, at closest (^a.b) and farthest (0 a,b) separation. In Fig. 7b we 
zoom in to see how the angles are defined^ 


fA,B = tan 


/ flA,B sin A/ -b B 
\ flp - Oa.B cos a/ 


( 22 ) 


and 


0a,b = tan 


aA,B sin A/ -b a.b 
flp -t- r^A.B cos a/ 


(23) 


The larger angle (is what corresponds to the limit of transitability 
and hence 6 does not appear in any further equations. 

All observers within the hatched band will eventually see tran¬ 
sitability. The probability of an observer being within the hatched 


(a) 




Figure 7. In (a) we show a side-on view of a circumbinary system with the 
planet shown at its two extreme values of /p. Both stars are plotted twice 
in different colours at closest and farthest separation from the planet. The 
angles subtended by the orbit on the celestial sphere are only shown for the 
primary star, to avoid clutter. The hatched region corresponds to observers 
who would see the planet in transitability at some point in time. In (b) we 
zoom in to see how d and f are defined. 


band in Fig. 7a is 


Fa.b = 


^ir/2 

Ji/2-A 


sin /biil^fbin 

nl2-hI-( 

= sin(A/ -b ^a,b) 

^A.B sin A/ -b /?a,b 


= sin A/ -b 


flp - Aa.b cos a/ 


(24) 


where because ^ is generally small we can apply the small angle 
approximation to remove the tan“* function'*. Equation 24 is the 
probability of transitability on the primary and/or secondary stars. 


^ The current equations are corrected from typos that appeared in the pub¬ 
lished version: in Eq. 22 there used to be a -b sign in the denominator and 
in Eq.23 there used to be a - sign in the denominator. We are sorry for the 
error. 


* Related to the eiror mentioned in footnote 3, Eq. 24 has been corrected 
from the published version which incorrectly contained an -b sign in the 
denominator. This error was purely typographical and all results shown in 
the paper were calculated using the correct equation. 
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for a binary of any orientation. The inclusion of f adds a period- 
dependency that is absent in Eq. 21. 

In Fig. 8a we demonstrate Eq. 24 on an example circumbinary 
system, comprised of a binary with Ma = IMq, = 

0.5Mo, = O.5I?0 and abin = 0.082 AU (Tbin = 7 days). The 

planet semi-major axis is varied from 0.24 AU to 2 AU. The three 
mutual inclinations are 0°, 5° and 10°. 

As A/ is increased the probability of transitability is increased 
significantly. For the misaligned cases, Pb > Pk- In Fig. 8h we 
zoom in on the coplanar case. As a comparison, we show the transit 
prohahility on a single star of radius Ra,b^ calculated using 

Ea,b = — . (25) 

For coplanar systems the probability of transitability reduces 
to 


Pks = . (26) 

Op - flA.B 

which comes from setting A/ = 0 in Eq. 24 and using a small angle 
approximation to remove the sin function. This equation matches 
Welsh et al. (2012), who derived an analytic estimate for the tran¬ 
sit probability^ under the assumption of Afl = 0 and static orhits, 
although it was duly noted that circumbinary orbits precess. Equa¬ 
tion 26 is close to the single star prohahility hut slightly higher, 
particularly at short periods. This is because the stars are brought 
closer to the planets by their orbital motion. 

The crucial difference to the single star case is that Eq. 25 de¬ 
creases towards zero for large semi-major axes, but for circumbi¬ 
nary systems the limit of Eq. 24 is 


lim Pa.b = sin A/, (27) 

flp—»oo 

which is equal to the first-order derivation in Eq. 21. This approx¬ 
imate probability is also applicable to systems with a large A/ be¬ 
cause the angle which encompasses the period-dependence, be¬ 
comes relatively small. 

In Fig. 9 we demonstrate how the probability of transitability 
varies with A/, using the same circumbinary system as in Fig. 8a, 
but fixing flp = 0.26 AU. In the bottom right of this figure we zoom 
in near A/ = 0°. The curve for the secondary star is seen to overtake 
that of the primary at around A/ = 3°. 


6 CONSEQUENCES FOR ECLIPSING BINARIES 

Eclipsing binaries are only a small fraction of the total binary popu¬ 
lation but the easiest binaries to detect photometrically. It was sug¬ 
gested by Borucki & Summers (1984) that they are favourable tar¬ 
gets for transit surveys because they positively bias the planetary 
orbit towards being aligned with the line of sight. In this section 
we derive the probability of transitability, Pa.b. under assumption 
that the binary is known to eclipse. In doing so we quantify what 
was first noted in Sect. 4 upon deriving Eq. 7: the criterion for tran¬ 
sitability is easiest to fulfil in the case of eclipsing binaries. The 
criterion for an eclipse is 


^ In fact, their derivation was for the probability of transitability, despite 
not using that name. 



Figure 9. The probability of transitability on stars A and B calculated using 
Eq. 24, and the approximation using Eq. 21. The bottom right image is 
zoomed into small mutual inclinations. In this plot the horizontal lines are 
the equivalent single star transit probabilities (Eq. 25). 


sin 


n 

2 


Rk + Rr 


■ 2.(xRr 


Obii 


(28) 


where a determines whether the criterion is for grazing eclipses 
(a = 0), full eclipses (a - 1) or anything in between. Since eclipses 
occur when /bin ~ ?r/2, we can apply the small angle approximation 
in Eq. 28 to deduce that the distribution of /yn for eclipsing binaries 
is uniform between 


where to simplify the equation we have defined 

S = RA + RB-2aRB. (30) 

Knowing that the binary eclipses with this uniform random 
distribution of /bin, the probability of transitability is the fraction 
of eclipsing binaries with /bin such that the inequality in Eq. 17 is 
satisfied. First, multiply Eq. 28 by Pa,b and insert it into Eq. 17, 
using /Ja.b from Eq. 19, to obtain 

A/ = Pa.b— - sin^‘ (f’A.B/^B.A —— - • (31) 

^bin \ ^bin ^p / 

The inequality from Eq. 17 has disappeared since we are calculat¬ 
ing Pa,b for a given A/. Use the small angle approximation and 
rearrange to obtain 


A/ = Pa.bS 
By solving for Pa,b we get 


rtbin 


■ Pb,a - 


1 \ Ra 


PA,i 


AI + 


R. 


A,B 




1 1 
-A^b,a — 

^bin 


if AI < A/iin, 


if AI > A/iin, 


(32) 


(33) 
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Figure 8. In (a) we show the probability of transitability on stars A and B (Eq. 24) as a function of dp, for three dilferent mutual inclinations. The horizontal 
dashed lines are calculated using the first-order approximate probability (Eq. 21). In (b) we zoom in on the Al = 0° case. As a compaiison, we show the 
equivalent single star probability in dashed lines (Eq. 25). 


where we define 

/I 1 \ ^AB 

A/i.m = <5- A<b,a- - —’ (34) 

\ C^bin CJp j Up 

in order to truncate Pa b at 1. As an example, for a binary with 
solar and half-solar masses and radii and abin = 0.3 AU orbited 
by a planet at Up = 1.0 AU, a = 0.5 and coplanar orbits, Eq. 33 
yields Pa = 0.33 and = 0.19. Coplanar orbits correspond to 
a minimum value of Pa.b. A slight increase in A/ to 0.5° raises 
these probabilities to 0.96 and 0.89, respectively. For A/ = 1° the 
probability on both stars is 1. The circumbinary geometry is evi¬ 
dently very favourable for transitability on eclipsing binaries. We 
note for reference that the mean mutual inclination in the transiting 
circumbinary planets found so far is 1.73° (see Table 1) and that 
the Solar System mutual inclination distribution roughly follows a 
1° Rayleigh profile relative to the invariant plane (e.g. Clemence & 
Brouwer 1955; Lissauer et al. 2011) 

For A/ > A/iin, two things occur: 1) transits are guaranteed on 
eclipsing binaries of any orientation and 2) transits become possible 
on non-eclipsing binaries. In Fig. 10 we plot A/nm as a function of 
flbin and Up, for a binary with solar and half-solar mass and radius 
and Qbiii between 0.007 and 0.2 AU, where the lower limit corre¬ 
sponds to a contact binary: Cbin = Ra+Pb- For the eclipse criterion 
we used a = 0.5. The white empty space on the left is the unstable 
region according to Eq. 4. 

The binary semi-major axis is the biggest factor in the calcu¬ 
lation. For the systems in Fig. 10 where circumbinary planets have 
been found so far (Obin > 0.08 AU), A/u^ is less than 3°. 

For closer binaries A/um rises sharply, reaching a maximum 
of 38° for a contact binary. Transits on very short-period eclipsing 
binaries are of course possible but Eq. 34 shows that not all such 
binaries can be transited unless there is significant misalignment. 
In Sect. 8.4 we apply this work to the Kepler discoveries so far. 

There are similarities between transits on eclipsing binaries 
and studies of multi-transiting systems orbiting single stars (e.g. 
Ragozzine & Holman 2010; Gillon et al. 2011). Geometrically, one 


is more likely to find planets transiting a single star where another 
transiting planet has already been found, compared to around a 
random star. This is because the mutual inclination distribution of 
multi-planet systems is not isotropic but weighted towards copla¬ 
narity*. This is analogous to a circumbinary system, if one con¬ 
siders the secondary star as the “inner planet”. There is, however, 
a fundamental difference between single and binary stars: planets 
orbiting single stars do so on effectively static orbits. 


7 CONNECTING TRANSITABILITY TO TRANSITS 

7.1 Does transitability guarantee transits? 

Transitability alone is not detectable via photometry, one requires 
an actual transit. A fundamental element of the definition of tran¬ 
sitability is that transits are possible but not guaranteed on any 
given passing of the binary orbit. By having a transit probability 
between 0 and 1 for each passing, it is intuitive to think that a tran¬ 
sit will eventually happen if observed continuously for a sufficiently 
long time. This conclusion was shared by Schneider (1994); Welsh 
et al. (2012); Kratter & Shannon (2013); Martin & Triaud (2014). 

We tested this hypothesis by numerically simulating circumbi¬ 
nary systems over 50 precession periods and looking for transits 
in cases where transitability occurred. This was first done for com¬ 
pletely random systems taken from the tests in Appendix A. The de¬ 
tails are provided in Appendix B. Less than 0.3% of systems man¬ 
aged to evade transit. All of these exceptional cases corresponded 
to the limit of transitability, where the planet only spends a very 
short time in transitability, and hence the chance of transiting on 


* Ragozzine & Holman (2010) found it resembles a Rayleigh distribution. 
The distribution of A/ in circumbinary systems is presently unknown, be¬ 
cause the detections so far have been highly biased towards coplanarity 
(Martin & Triaud 2014). 





















8 Martin and Triaud 



Figure 10. 3D histogram of the minimum mutuai inciination needed in de¬ 
grees to guarantee transits on EBs of any orientation, at different binary and 
pianet semi-major axes (Eq. 34). 


any given orbit is small. Transits are expected to occur eventually, 
but after a time longer than what was simulated. 

Second, we constructed systems with 4:1 and 5:1 period com- 
mensurabilities, specifically designed to make planets permanently 
evade transit. The evasion percentage increased slightly but re¬ 
mained less than 1%. It is probable that this value would eventually 
drop to zero, but longer simulations would be required. The planets 
inevitably transit because exact period commensurabilities are not 
sustainable, owing to perturbations from the binary on the plane¬ 
tary orbit (Sect. 8.5.2). Aside for HD 202206 which has a period 
ratio near 5:1 (Correia et al. 2005), period commensurabilities have 
not been observed. This system may not be representative, since it 
straddles the border between a circumbinary and a two-planet sys¬ 
tem; the secondary “star” has a minimum mass of 15Mj„p. It has 
also been theorised by Kley & Haghighipour (2014) that circumbi¬ 
nary planets should form between integer period ratios, not in them. 

Whilst not an exhaustive proof, our tests indicate that in the 
vast majority of cases, transitability indeed leads to transit, albeit at 
an unspecified point in time. 

7.2 Transits over time 

The probability of transitability is equivalent to the probability of 
transit, granted the observer has infinite time. Unfortunately, due 
to limitations in technology, funding and human life-expectancy, 
one must strive to capture transits within a finite time. We calcu¬ 
late some example observation times needed in order to observe a 
transit. 

In the case of a single star, after continuous observations of a 
time equal to Tp, the planet either will or will not have transited. 
This is not the case for circumbinary planets, for two reasons: 

• The planet may currently be outside of transitability, but will 
precess into transitability at a later time. 

• The planet may currently be inside transitability, but the con¬ 
junction required for a transit has not yet occurred. 

The fraction of circumbinary planets transiting therefore in¬ 
creases with time, up to a value specified by Eq. 24. It is important 


cn 



Figure 11. The percentage of systems seen transiting as a function of time 
from numerical simulations of 10,000 circumbinary systems, and the ana¬ 
lytic predictions (Eq. 24). Systems are counted as transiting after the detec¬ 
tion of a single transit. The zoomed figure in the bottom right comer shows 
the percentage of transiting systems over kep/er-like observing times. The 
black vertical line denotes the precession period. 


to know how long an observer must wait to see a transit. It is a 
strong function of the precession period, since that determines how 
spaced apart the regions of transitability are. 

An analytic calculation of the time-dependent transit proba¬ 
bility is outside the scope of this paper, and has been previously 
labelled impossible (Schneider & Chevreton 1990). We instead use 
numerical N-body simulations. 

In Fig. 11 we demonstrate the percentage of systems seen tran¬ 
siting stars A and/or B as a function of time, using 10,000 simulated 
circumbinary systems. The primary and secondary stars are solar 
and half-solar in mass and radius, Tbin = 7 days, Tp = 40 days and 
A/ = 10. Over time, the percentage of transiting systems reaches 
the value predicted by Eq. 24, in agreement with the conclusions 
of Sect. 7.1. Most of the transiting systems have done so within a 
single precession period (here ~ seven years). 

As an extended test, we took the systems found transiting in in 
Sect. 7.1 and calculated the time taken for primary and secondary 
transits to occur. The results are provided in Table Cl. Whilst a 
larger mutual inclination leads to more planets transiting, the me¬ 
dian wait time is increased. Generally, a significant number of sys¬ 
tems are found transiting within A'ep/er-like mission times. 


8 DISCUSSION AND APPLICATIONS 

8.1 The circumbinary planets discovered so far 

Our first application is to calculate the transitability probabilities 
for the Kepler discoveries so far, assuming of course that we do 
not have a priori knowledge of transits and eclipses^. In Table. 1 
we calculate the probability of transitability on binaries of any ori¬ 
entation (Eq. 24) and on eclipsing binaries (Eq. 33), where for the 
latter we used a = 0.5 to define eclipses. The equivalent single star 

^ Otherwise you would have a boring table full of 100%’s. 
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Table 1. Probabilities of transit for the circumbinary planets detected so far by Kepler. 


Name 

Ma 

(Me) 

Mb 

(Me) 

Ra 

(Re) 

Rb 

(Re) 

^bin 

(AU) 

(AU) 

A/ 

(deg) 

Ra.b 

A 

% (all) 

B 

Ra.b 

A 

% (EBs) 

B 

Ra.b 

A 

% (single) 

B 

Kepler-16 

0.69 

0.20 

0.65 

0.23 

0.22 

0.71 

0.31 

1.04 

0.91 

75.5 

66.1 

0.42 

0.15 

Kepler-34 

1.05 

1.02 

1.16 

0.19 

0.23 

1.09 

1.86 

4.18 

4.16 

100 

100 

0.50 

0.47 

Kepler-35 

0.89 

0.81 

1.03 

0.79 

0.18 

0.60 

1.07 

3.11 

2.94 

100 

100 

0.78 

0.61 

Kepler-38 

0.95 

0.27 

1.78 

0.27 

0.15 

0.46 

0.18 

2.28 

0.79 

41.2 

14.4 

1.80 

0.27 

Kepler-47b 

1.04 

0.46 

0.84 

0.36 

0.08 

0.30 

0.27 

1.93 

1.26 

39.5 

25.8 

1.30 

0.56 

Kepler-47c 

1.04 

0.46 

0.84 

0.36 

0.08 

0.99 

1.16 

2.48 

2.32 

50.8 

47.6 

0.39 

0.17 

Kepler-64 

1.50 

0.40 

1.75 

0.42 

0.18 

0.65 

2.81 

6.54 

6.66 

100 

100 

1.25 

0.30 

Kepler-413 

0.82 

0.52 

0.78 

0.48 

0.10 

0.36 

4.02 

8.98 

9.18 

100 

100 

1.01 

0.62 

KIC9632895 

0.93 

0.19 

0.83 

0.21 

0.18 

0.93 

2.30 

4.68 

5.09 

100 

100 

0.49 

0.12 


Refs: Doyle et al. (2011); Welsh et al. (2012); Orosz et al. (2012a, b); Schwamb et al. (2013); Rostov et al. (2013, 2014) 
Welsh et al. (2014) 

Note: Kepler-47d is excluded because it has not yet been published and lacks a value for A/. 


probability was calculated using Eq. 25. In the table we include all 
necessary variables for the calculations. In more than half of the 
cases, transits are guaranteed on eclipsing binaries of any orienta¬ 
tion. 

8.2 Multi-planet circumbinary systems 

Only one multi-planet circumbinary system has been discovered so 
far (Kepler-47, Orosz et al. 2012b). Kratter & Shannon (2013) con¬ 
sidered an eclipsing binary with a known transiting planet, and cal¬ 
culated the likelihood of a second planet being seen transiting. They 
derived an analytic probability for whether or not the binary and 
planet orbits would overlap on the sky, under the assumption that 
the binary is perfectly edge-on (/bin = ^/2). In fact, what they cal¬ 
culated was the probability of transitability. Their derivation does 
not include precession, and consequently underestimates the prob¬ 
ability. 

Based on the work in Sect. 6, any additional planets with A/ 
greater than the first transiting planet are guaranteed to enter tran¬ 
sitability at some point. 

8.3 Kepler’s eclipsing binary catalog 

The Kepler telescope, with its four years of continuous observa¬ 
tions and exquisite precision, has provided the most comprehen¬ 
sive catalog of EBs to date (Slawson et al. 2011). We used the 
online beta version of this catalog* to test our transitability crite¬ 
rion on hypothetical orbiting planets. From the catalog we obtained 
M[, M 2 ,Ri,R 2 and flbin. which were derived from stellar temper¬ 
atures calculated in Armstrong et al. (2013) using a method ex¬ 
plained in Armstrong et al. (2014). Only systems with a morphol¬ 
ogy parameter less than 0.5 were used, corresponding to detached 
EBs (see Matijevic et al. 2012 for details). The binary inclination 
was randomised between the bounds defined in Eq. 29^. 

The remaining quantities needed for Eq. 17 are Al and flp. 
Given the distribution of circumbinary planets is presently poorly 
known and subject to strong biases, we considered a wide range of 
potential values. In Fig. 12 we calculated the percentage of Kepler 

* http://keplerebs.villanova.edu/ maintained by Andres Prsa et al. 

^ It is possible to obtain a tme value of /yn, however the only published 
version is in the now-outdated catalog of Slawson et al. (201 1), and contains 
errors. 


EBs on which there would be transitability by a putative planet, 
with Al between 0° and 10° and Qp/abin = 3,10,20. The results are 
only shown for the primary star, since the plot for the secondary 
star is indistinguishable. 

Consistent with earlier sections, the biggest factor is the mu¬ 
tual inclination, with a higher Al leading to a greater chance of 
transitability. Transitability is favoured for smaller values of dp, but 
this dependency diminishes at larger mutual inclinations. 

The Kepler EB catalog would benefit from extended photo¬ 
metric observations by the future PLATO telescope (Rauer et al. 
2014), in order to find new circumbinary planets that have moved 
into transitability during the ~ eight years between missions. This 
may include additional planets in known circumbinary systems. 


8.4 On the dearth of planets around short-period binaries 

An observed trend has been the lack of circumbinary planets around 
the closest binaries; the shortest-period binary hosting a planet is 
Kepler-47 with Tbin = 7.4 days (Orosz et al. 2012b). This is de¬ 
spite the median*® period of the EB catalog being 2.8 days. This 
raises various questions about the ability to form planets in such an 
environment, particularly in the presence of tertiary stellar compan¬ 
ion, as is often the case for very tight binaries according to theory 
(Mazeh & Shaham 1979; Fabrycky & Tremaine 2007) and obser¬ 
vations (Tokovinin et al. 2006). 

The reason why EBs are preferentially found at short peri¬ 
ods despite a smaller natural occurrence (Tokovinin et al. 2006) 
is because there is a greater range of /bin that allow for an eclipse 
(Eq. 29). When the EB is highly inclined, however, the planet itself 
needs a greater misalignment in order for transitability to occur. 
We demonstrate this in Fig. 13, where we calculate the minimum 
Al needed to see transitability on the primary star of each EB in the 
Kepler catalog, taking ap/flbin = 3.5. As a reference, we show the 
mean and maximum mutual inclinations from the Kepler discover¬ 
ies so far, although these are biased towards being small (Martin & 
Triaud 2014). 

For Tbin < 7.4 days, a misalignment of 1.73° results in 60% of 
planets missing transitability. For those systems misaligned enough 
for transitability on the shortest period binaries, there should be 


'® In the published version of the paper we accidentally wrote that this was 
the “mean” period. 
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Figure 12. The percentage of EBs found by Kepler on which there would Figure 13. The minimum mutual inclination needed for transitability on 

be transitability by a putative planet with different values of A/ and Op/abin. each of the EBs found by Kepler. 


transits within the Kepler time series, for two reasons: 1) the pre¬ 
cession period is only a couple of years long", so the planet and 
binary orbits would have intersected at least once during the Ke¬ 
pler mission and 2) these are very tight systems, so we are likely to 
have observed one or probably more transits whilst in transitability. 

The dearth of planets may also be explained by stellar noise 
in the light curves. Binaries with periods this short are expected to 
be tidally locked, which leads to faster rotation and increased star 
spots, which may inhibit detections. 

The current null detection likely remains significant for plan¬ 
ets that are misaligned by at least a few degrees, but there may re¬ 
main some coplanar ones that are undetectable by Kepler. Transits 
on contact binaries require an even higher level of mutual inclina¬ 
tion. Discoveries around contact binaries may also be hindered by 
Kepler’s 30-minute cadence, which is potentially too long to ade¬ 
quately sample its orbit in the search for transits. 

8.5 Limitations 

8.5.1 Eccentric Systems 

The addition of eccentricity, to both the binary and planet orbits, in¬ 
troduces two complexities. First, the geometry is complicated since 
we lose circular symmetry, and there are two additional angles to 
consider: Wbin and ojp (the arguments of periapse). Furthermore, the 
orbital dynamics cause oip to be time-dependent, further complicat¬ 
ing the situation. 

Second, the precession cycle is more complex when the binary 
is eccentric. In Fig. 14 we demonstrate the precession of the same 
system as in Fig. 3, but with gbin = 0.5. The mutual inclination is 
no longer constant. There are two islands of libration, centred on 
AD = 0 and A/ = n/2 (red) and A/ = -n/l (magenta), within 
which AD does not circulate through 0 to 2n. We therefore lose 
two of the assumptions made in Sect. 4. 

We ran numerical simulations to test the ability of our criterion 
in Eq. 17 to predict transitability, in the case of eccentric systems. 

" If we assume the observed over density of planets at Tp ~ STbin extends 
to very close binaries. 



Figure 14. The same circumbinary systems as in Fig. 3 but with ebin = 0-5. 


The simulations show that the error is only of order ~ 2% (Ap¬ 
pendix D). Furthermore, the results suggest that eccentricity actu¬ 
ally makes transitability more likely, although a more detailed study 
is needed to confirm this. 

8.5.2 Additional dynamical effects 

A circumbinary system is a three-body problem and hence not solv¬ 
able analytically. Our analytic treatment of it as a pair of Keplerians 
plus orbital precession encompasses the majority of the physics, but 
neglects some smaller amplitude effects. 

The semi-major axis and period, which we assumed to be con¬ 
stant, experience slight variations over time because of perturba¬ 
tions from the binary. There are also small variations in the ec- 
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Figure 15. Osculating orbital elements for an example circumbinary sys¬ 
tem, viewed over a timespan slightly greater than 2/’p. 


centricity. In fact, even a system with an initially circular planetary 
orbit will obtain some eccentricity over time, and hence the compli¬ 
cations of Sect. 8.5.1 are to a certain extent unavoidable. In Fig. 15 
we present an example of the variation in planet period (top), semi¬ 
major axis (middle) and eccentricity (bottom) for a binary with so¬ 
lar and half-solar mass and radius stars, Tbin = 10 days and a planet 
with an initial period of Tp = 52 days and an initially circular orbit. 
There is a short-term variation with period l/TTbin and a longer- 
term modulation with period Tp. A comprehensive analysis of the 
orbital dynamics of circumbinary systems can be found in Leung 
&Hoi Lee (2013). 

We have also neglected any effects that may be imposed by 
a tertiary star. Tertiary stars are very commonly found around bi¬ 
naries (Tokovinin et al. 2006) and, if close enough, may affect the 
precession and observability of any circumbinary planets. It is also 
possible that tertiary stars with sufficient influence would also have 
hindered planet formation; this case would consequently be rare. 

Additional dynamical effects may also arise in the presence of 
multi-planet circumbinary systems, although it remains to be seen 
if these would have a noticeable impact in comparison with the 
perturbations from the binary. 


8.5.3 The detectability of transits 

Whilst we have shown that transits are more likely at a higher mu¬ 
tual inclination, we have not considered their detectability. Even in 
the simplest case of a coplanar planet, there are significant varia¬ 
tions in transit timing (Armstrong et al. 2013) and duration (Liu et 
al. 2014), making detections more difficult than in the single star 
case. For highly misaligned systems, the transits will be aperiodic, 
and perhaps even singular (Martin & Triaud 2014). Also, because 
of the precession, the transit signature may disappear for a while. 
It was not until the discovery of the ninth transiting circumbinary 
planet - the 4.1° misaligned Kepler-413 by Kostov et al. (2014) - 
that the effects of time-dependent transitability were readily appar¬ 
ent. This does demonstrate, however, that detection techniques are 
improving in order to discover these more complicated systems. 

Another unsolved question is the predictability of transits. 
Martin & Triaud (2014) showed that the sequence of primary and 
secondary transits observed for misaligned systems is highly sensi¬ 
tive to the input orbital parameters, e.g. altering the starting planet 
inclination by as little as 1° may completely change the transit num¬ 
ber and timing. The small variations in orbital elements described 
in Sect. 8.5.2 will have to be accounted for, although the regular¬ 
ity of the variation seen in Fig. 15 is promising. Predictability is 
not a concern for blind, continuous surveys like Kepler, TESS and 
PLATO. It is crucial, however, for any targeted follow-up observa¬ 
tions. 


9 CONCLUSION 

The geometry and orbital dynamics unique to circumbinary planets 
make them very likely to exhibit transitability, and hence transit at 
some time. Transits provide exoplanetary science with a wealth of 
information. This probability can be increased to one when observ¬ 
ing eclipsing binaries. Furthermore, there is a relatively weak de¬ 
pendence on period. This allows an extension of transit surveys to 
the outer regions of exoplanetary systems, with applications includ¬ 
ing including the atmospheric characterisation of cold exoplanets. 
The high transit probability bodes well for complementary obser¬ 
vations of circumbinary planets found with other techniques, such 
as ETVs (Borkovits et al. 201 1), radial velocimetry (Konacki et al. 
2009) and astrometry (Salmann et al. 2014). For these reasons, cir¬ 
cumbinary planets can be seen as more than just exotic examples 
of nature’s diversity but as practical tools in astronomy. 
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APPENDIX A: TESTING THE ANALYTIC CRITERION 

The test consisted of creating 10,000 random circumbinary sys¬ 
tems, integrating them over an entire precession period*^, and 
checking if the planet and stellar orbits overlapped at any time. The 
orbital parameters were drawn from a uniform distribution between 
the maximum and minimum bounds listed below. 

• For the binary: Ma between 0.5 and 1.5 Mq; Mb calculated 
using a random mass ratio q drawn between 0.2 and 1.0; Ra and 
calculated using a mass-radius relation of Kippenhahn & Weigert 
(1994); Tbin between 5 and 50 days; fibin = 0 and the starting orbital 
phase between 0 and In. 

• For the planet: Mp = 0; Rp = 0; Tp/Tbin between 4 and 10; 
gp = 0 and the starting orbital phase between 0 and In. 

• Two uniform mutual inclination distributions were tested: 0° 
to 5° (test 1) and 0° to 30° (test 2). We drew AO between 0 and 
In. The orientation of the circumbinary system on the sky was ran¬ 
domised using a uniform 3D rotation algorithm by Arvo (1992), 
which creates a uniform distribution of cos /bbi. 

These orbital parameters are not completely arbitrary, since 
they encompass all of the circumbinary systems discovered to date. 
The planet period was chosen with respect to the binary period so 
that the stability limit was respected (Eq. 4). 

The test results are shown in Table AL The result for each 
system was put into one of four categories: 

• Cat. 1: Eq. 17 predicted transitability and the numerical simu¬ 
lation matched this. 

• Cat. 2: Eq. 17 predicted no transitability and the numerical 
simulation matched this. 

The integration time was LlTprec, where Tprec was calculated using 
Farago & Laskar (2010), and the factor of 1.1 allows for any errors in their 
formula. 
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• Cat. 3: Eq. 17 predicted transitability but the numerical simu¬ 
lation did not show transitability 

• Cat. 4: Eq. 17 predicted no transitability but the numerical 
simulation showed transitability 

The last two categories correspond to errors in the analytic 
formula. 


APPENDIX B: TESTING THE CONNECTION BETWEEN 
TRANSITABILITY AND TRANSITS 

As a starting point, we took the circumbinary systems from Ap¬ 
pendix. A that were found in transitability (test 1: A/^ax = 5° 
and test 2: A/max = 30°). The small number of error cases in Ap¬ 
pendix. A were avoided. For each system we integrated over 50 
precession periods to look for primary and secondary transits. 

Next, we tested whether planets on periods commensurate 
with that of the binary would manage to permanently evade transit. 
We created 10,000 systems with the same parameters as in Ap¬ 
pendix A, but with Tp/Tbin fixed at 4 and at 5. This led to four 
additional tests: 3) A/^ax = 5° and Tp/Tbin = 4; 4) A/^ax = 5° and 
Tp/Tbin = 5; 5) A4ax = 30° and Tp/Tbia = 4; and 6) A/„ax = 30° 
and Tp/Tbin = 5. 

In Table B1 we show the number of systems in transitability 
and the number of systems evading transit during 50 precession 
periods, for the six tests outlined above. 


APPENDIX C: MEDIAN TRANSIT WAIT TIMES 

For each of the transiting circumbinary systems in Appendix B we 
calculated the time of first transit on each star. In Table. Cl we show 
the median time taken. We also show this time as a fraction of the 
precession period. In the last two columns we show the percent¬ 
age that are seen transiting at least once within four years, which 
comparable to a Kepler-Xtka mission. 

The purposes of these numbers is to provide a rough estimate 
of the time needed for transits to occur, and to motivate the fact that 
they can occur within realistic timeframes. 


APPENDIX D: ECCENTRIC SYSTEMS 

We ran a single set of numerical simulations, similar to in Ap¬ 
pendix A, to test how accurate the transitability criterion was for 
eccentric systems. All of the parameters used were the same as in 
Appendix A, except Cbin and Cp were varied uniformly between 0 
and 0.5. The mutual inclination maximum was also increased to 
60°, so we could probe the islands of libration seen in Fig. 14. 
Adding eccentricity can make a system unstable, so we removed 
any such cases according to either the stability criterion in Holman 
& Wiegert (1999). Since Holman & Wiegert (1999) includes ebin 
but not Cp, we used the periapsis distance ap(l-ep) as a approximate 
means of including the effects of planet eccentricity. Additionally, 
we monitored and removed planets that were ejected during the nu¬ 
merical simulations. The results are shown in Table. Dl. 

In Fig. Dl we plot a histogram of A/ for the 105 systems in 
Cat. 3, for both A and B stars, where transitability was predicted 
but did not occur. There is a large jump at roughly 50°. This corre¬ 
sponds to the islands of libration (Fig. 14). Such high mutual incli¬ 
nations are predicted to yield transitability according to Eq. 17, but 



A I (deg) 


Figure Dl. Histogram of the 105 systems in Cat. 3 for which transitability 
was predicted on either star in an eccentric system but did not occur. 


since Afl no longer passes through all angles, some unlucky aliens 
will never see transitability. 
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Table Al. Testing the transitability criterion 


test A/max (deg) Cat. 1 Cat. 2 Cat. 3 Cat. 4 error % 


A B A BABABA B 

1 5 618 643 9381 9352 1 2 0 3 0.01 0.03 

2 30 2884 3105 7106 6884 9 4 1 4 0.06 0.06 


Table Bl. Testing if transitability leads to transits 


test A/max (deg) Tp/Tbin no. in transitability no. evading transit evasion % 





A 

B 

A 

B 

A 

B 

1 

5 

4- 10 

618 

641 

0 

2 

0 

0.31 

2 

30 

4- 10 

2878 

3096 

7 

9 

0.24 

0.29 

3 

5 

4 

621 

663 

3 

4 

0.48 

0.60 

4 

5 

5 

593 

662 

2 

6 

0.34 

0.91 

5 

30 

4 

3044 

3370 

15 

30 

0.49 

0.89 

6 

30 

5 

2956 

3240 

11 

14 

0.37 

0.43 


Table Cl. The time taken for primary and secondary transits to occur 


test A/max (deg) Tp/T^in median time (yr) scaled median time percentage < 4 yr 





A 

B 

A 

B 

A 

B 

1 

5 

4- 10 

6.6 

8.4 

0.24 

0.31 

37 

33 

2 

30 

4- 10 

18.5 

28.0 

0.50 

0.74 

19 

16 

3 

5 

4 

2.4 

3.4 

0.24 

0.35 

68 

55 

4 

5 

5 

3.4 

4.5 

0.22 

0.30 

55 

46 

5 

30 

4 

6.9 

11.7 

0.63 

1.03 

37 

27 

6 

30 

5 

10.2 

16.5 

0.57 

0.85 

29 

21 


Table Dl. Testing the transitability criterion on eccentric systems 


A/max (deg) Cat. 1 Cat. 2 Cat. 3 Cat. 4 Error % 


A B A BABABAB 
60 3185 3376 2677 2493 96 65 146 170 1.78 1.72 



